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By using analytical solution of a tight-binding model for armchair nanoribbons, it is con- 
firmed that the solution represents the standing wave formed by intervalley scattering and 
that pseudospin is invariant under the scattering. The phase space of armchair nanoribbon 
which includes a single Dirac singularity is specified. By examining the effects of boundary 
perturbations on the wave function, we suggest that the existance of a strong boundary po- 
tential is inconsistent with the observation in a recent scanning tunneling microscopy. Some of 
the possible electron-density superstructure patterns near a step armchair edge located on top 
of graphite are presented, ft is demonstrated that a selection rule for the G band in Raman 
spectroscopy can be most easily reproduced with the analytical solution. 
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1. Introduction 

Graphene nanoribbons are particularly useful in study- 
ing the properties which are hidden in carbon nanotubes 
by the periodic boundary condition. 1 ) The boundary of 
a nanoribbon consists of two symmetrical edge struc- 
tures, armchair and zigzag edges. A transmission electron 
microscope image illustrates that the edge of a sample 
prepared by micromechanical cleavage, 2 ) is a mixture of 
armchair and zigzag edges. 3 ) Some attempts have been 
made to prepare a regular zigzag or armchair edge over 
a large length more than 20 mn. 4 ' 

It is known that the behavior of the electrons near 
zigzag edge differs greatly from that near armchair edge 
in many respects. For example, armchair edge gives rise 
to intervalley scattering, while zigzag edge produces in- 
travalley scattering. This originates from the fact that 
the graphene Brillouin zone in the /c-space is given by 
rotating the hexagonal unit cell in real space by 90°. 5 ) 
On the other hand, zigzag edge changes the orientation of 
pseudospin, while armchair edge does not. This behavior 
of pseudospin results from that only A or B (both A and 
B) sublattice appears at zigzag (armchair) edge. These 
differences between zigzag and armchair edges cause ob- 
servable effects to appear in Raman spectra. Because the 
D band in Raman spectra is relevant to intervalley scat- 
tering, it follows that the D band intensity is enhanced 
by armchair edge. 6 ) The behavior of pseudospin is essen- 
tial to a selection rule of the G band for the graphene 
edge. 7 ) 

Besides the differences in the behavior of the scattering 
and pseudospin, armchair and zigzag edges are distinct 
with respect to the number of Dirac singularities. In the 
case of nanotube, there are two Dirac singularities corre- 
sponding to the K and K' points at the corners of the first 
Brillouin zone. By means of an effective-mass theory, we 
discussed in a previous paper that there is a single Dirac 
singularity in the phase space of armchair nanoribbon, 



while there is no Dirac singularity in the phase space of 
the standing wave of zigzag nanoribbon. 5 ) The existence 
of a single Dirac singularity in armchair nanoribbon is 
of particular interest in transport property because the 
existence of the Dirac singularity can suppress backward 
scattering. The appearance of the edge states 8,9 ) near 
zigzag edge is a consequence of the absence of the Dirac 
singularity for zigzag nanoribbon. 

In this paper, by constructing analytical solution of 
a tight-binding model for armchair nanoribbons, we ac- 
count for the features of armchair edge, such as interval- 
ley scattering, the behavior of pseudospin, and the sin- 
gle Dirac singularity in the phase space. Since the tight- 
binding model has an atomic resolution, we can apply the 
solution to some recent experimental results obtained by 
scanning tunneling microscopy. 10,11 ) 

Here, we would like to mention the previously pub- 
lished literature on analytical solution for the wave func- 
tion in armchair nanoribbon. Compernolle et al. 12 ^ con- 
structed the wave functions in nanoribbons (achiral nan- 
otubes with edges) using the transfer matrix method. 
They applied the result to obtain the transmission co- 
efficient. Zheng et a/. 13 ) conducted analytical study of 
electronic state in armchair nanoribbons. They showed 
that all armchair nanoribbons acquire nonzero energy 
gap due to the variation of hopping integral near the 
edges. The above-mentioned authors employed a sim- 
ple nearest-neighbor tight-binding Hamiltonian that was 
shown to be a good approximation for describing the 
electrons in nanotubes and graphene systems. In this pa- 
per, we use the same Hamiltonian to examine armchair 
nanoribbons. 

This paper is organized as follows. In §2 we solve a sim- 
ple tight-binding model for armchair nanoribbons, and 
obtain the energy dispersion eq. (13) and the wave func- 
tion eq. (24) . The characteristic features of the wave func- 
tion are pointed out. By using eq. (24) we will show in 
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§3.1 that a feature of the Dirac singularity for armchair 
edge has been observed in recent scanning tunneling mi- 
croscopy experiment by Yang et al., 10 ^ and in §3.2 that 
the electron-phonon matrix element for the G band in 
Raman spectra which we have derived numerically with 
a tight-binding model in a previous paper 14 ) can be de- 
rived analytically. Discussion and summary are given in 
§4- 

2. Construction of Wave Function 

Let 4> J A and 4> B be the probability amplitudes of the A 
and B-atoms at the box in Jth line shown in Fig. 1(a), 
then the Schrodingcr equation for the electron in an arm- 
chair nanoribbon is written as the recurrence equation, 



+ 





-,+ikb 



+ 



tkb 



fr 1 



(i) 



Here, e is the energy eigenvalue in units of the hopping 
integral (7 = 3 eV) , k is the wave vector along the edge, 
and b (= 21) denotes the unit length along the axis of 
armchair nanoribbon [see Fig. 1(a)]. In obtaining cq. (1) 
we have used the Bloch theorem, by which the amplitude 
of the B-atom in J — 1th line that is located nearest to 
the A-atom at the box in Jth line is given by e~ lkh 4> J ^~ 1 . 
Note that J in cq. (1) takes from J = 2 to N — 1. For 
J = 1 and J = N, the recurrence equations are given by 

'° ^ W + ( r 

1 UU + \e+ ika 0, 
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When we assume that eq. (1) is satisfied for J = 
1, • • • ,N, these equations at J = 1 and J = N can be 
included as the boundary condition; 

0. x 1 .!) 0. (4) 



In the following, we solve the recurrence equation of 
eq. (1) with the boundary condition of cq. (4). 
First, let us introduce the following matrices, 



K = 



T = 



where z = e +lkb , z* = e 
eq. (1) is rewritten as 

T(j) J+1 + K(p J + T~ 1 (j) J - 





-tkb 



T -i = 



(5) 



, and kb G [— x, it). Then 



X =0, (J =!,•••, N) (6) 



where we have defined <j> J = t (cf>j i ,(j>' B ). By multiply- 
ing eq. (6) by T~ x from the left-hand side, we get 
4> J+1 + T- l K(f) J + (T- 1 ) 2 ^'- 1 = 0. Because the matrix 
T^ 1 satisfies (T -1 ) 2 = z*a , where a is 2 x 2 identity 
matrix, the recurrence equation can be diagonalized by 




(b) 




Fig. 1. (a) The structure of an armchair graphenc nanoribbon. 
The index J takes from 1 to AT. The lattice constant is a 
(= \/3a cc ). The unit length and the wave vector along the axis 
are denoted by b (= 3a cc ) and k, respectively. Carbon atoms are 
divided into A (•) and B (o) atoms, (b) The electron-density pat- 
tern of the Dirac singularity state in a metallic armchair nanorib- 
bon with N = 8. See also the STM topography shown in Fig. 3 
of Rcf. 10. 



some matrix U as 

tp J+1 + U~ 1 T~ 1 KUip J 



j-i 



= 0, 



(7) 



where cp J = U 1 (f> J . It is straightforward to show that 
the eigenvalues of the matrix T~ 1 K are given by 



A± 



l + z± v/(-l + z) 2 + 4ze 2 



2z 



and the corresponding eigen spinors are 



l-z±^/(-l+z) 2 +4ze 2 
1 



(8) 
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The eigenvalue A + (or A_) is a multivalued function since 
y/z has two branches: 



w + = +e 



ikl 



W- 
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(10) 



The eigenvalue A+ with the w+ (w-) branch is equiva- 
lent to the eigenvalue A_ with the W- (w + ) branch. As a 
result, we may consider the eigenvalue A + only. Here we 
consider the w + branch. We will show later that the W- 
branch gives the same wave function as the w + branch, 
and that these w + and W- branches correspond to dif- 
ferent states in the case of carbon nanotubes. 

Next, we solve the characteristic equation with respect 
to eq. (7) as 



x 2 + X + x + — 

z 

<^ (x — w^e 



= 
- i6 ) {? 



0, 



(11) 



where the parameter 9 e (0,7r) is defined by 2 cos 9 = 
—w + \ + or 

-7 1 + w + + \/z + z* - 2 + 4e 2 



2 cos 9 = 



(12) 



The solution to the characteristic equation is given by 
x = w^ 1 e + * e ' or x = w^e^ 16 . Since w~ [ 1 +w + = 2 cos(kl) 
and z + z* = 2cos(2fcZ), eq. (12) leads to the energy 
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dispersion relation, 

e 2 = 1 + 4cos 2 (9 + 4cos6lcosO0- (13) 

By employing the standard method for solving the recur- 
rence equation, general Jth term is written by the first 
and second terms, and as 

j = (J - 2) sin(J-l)fl _ (J _ 1)S in(J-2)fl 
sin# sin# 

(14) 

To this point, we have not yet used the boundary con- 
dition. The remaining procedure for obtaining the solu- 
tion is applying the boundary condition to eq. (14). The 
boundary condition of armchair edge is given by <f> = 0. 
Thus, eq. (7) with J = 1 gives <p 2 = -T^Kcj) 1 , so that 
we have 



of eq. (24) in terms of exponential functions as 



±2 x il 2cos# , 

By putting this into eq. (14), we obtain 
j _(j_i) sin J6 ! 

The eigen spinor <p\ may be represented in a simpler 
fashion as 



(15) 



(16) 
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where the variable O is a function of k and 9, and s takes 
+1 (—1) for a conduction (valence) band state [e = s|e|]. 
Therefore, the wave function is written as 



Ce -%ki(j-\) 



i JO 
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where C is the normalization constant. 

Several characteristic features of armchair nanorib- 
bons can be derived from the wave function eq. (18) and 
energy dispersion relation eq. (13). They are listed as 
follows. 

(I) Armchair nanoribbons are metallic for the case that 
N = 3i — 1, 15 ~ 17 ) where i is integer. From the energy dis- 
persion relation of eq. (13) one can see that the state 
with e = corresponds to the point (0, 27r/3) in the pa- 
rameter space (k,6). The parameter 9 is quantized as 
9 n = nn/(N + 1) [n = 1, • • • , N] by imposing the bound- 
ary condition on the wave function: (f>^ +1 (s, k, 9) = 0. 
The quantized parameter 9 n takes 27r/3 for the case that 
N = 3i — 1 and n = 2i. 

(II) The angle Q(k,9) in eq. (18) is singular at the 
point (0,27r/3). By putting 2tt/3 + 59) into 1 + 
2cos9/w+ of eq. (17), we h ave |e|e" ie = (-s/350 + iSkl) 
and |e| = yj?>59 2 + (Ski) 2 . The parameter 9 changes 
from to 27r when (5k, 59) moves along a path enclosing 
the singularity point. This behavior of results in a non- 
trivial Berry' phase which is a necessary condition for an 
armchair nanoribbon to exhibit a ballistic transport. 5 ' 18 -* 

(III) The wave function can be understood as the 
standing wave formed by intervalley scattering. This fea- 
ture can be made more transparent by rewriting sin( J9) 



f>+(s, k, 9) oc e 



-iklj 



-i@ 



—s 



— e 



-ije 



—s 



(19) 

where the parameter 9 e (0,7r) is a positive value. The 
exponential function e +lje (e~ lje ) indicates the propa- 
gating wave in the direction of increasing (decreasing) 
J. Hence, for the armchair edge at J = 1, the second 
term represents the incident wave to the edge and the 
first term represents the reflected wave. The momentum 
transfer through this scattering process, —9 — > 9, is given 
by 29. Since 9 = 2ir/3 (= 9p) is coincident with kpa/2, 
where fcp is the wave vector at the K point, we can see 
that the change in the wave vector at the armchair edge 
is given by 29-p = fcpa. This change of the wave vector 
corresponds to the intervalley scattering. 

(IV) The spinors of the incident and edge reflected 
waves in eq. (19) are equal. Hence, the pseudospin does 
not change under the intervalley scattering at armchair 
edge. This feature originates from the fact that both the 
A and B-atoms appear equivalently at the armchair edge. 
It is also clear that the probability densities at the A and 
B-atoms in Jth line are equal and the electron-density 
pattern depends only on the distance from the armchair 
edge. 

(V) There is a single Dirac singularity in the system. 
To make this point clear, we rewrite the standing wave 
eq. (19) as 



4>'+ = IpK - IpK' 



(20) 



where (tpK') corresponds to a Bloch state near the 
K (K 7 ) point. The minus sign in front of V'K' shows that 
this standing wave is an anti-symmetric combination of 
the two Bloch states. If the two states were independent 
of each other, we may expect that the symmetric combi- 
nation + V'K' is also a solution. However, such sym- 
metric solution should be excluded because V'K + V'K' is 
proportional to cos(J9) and does not satisfy the armchair 
boundary condition. Since the armchair boundary con- 
dition selects only the antisymmetric solution, the two 
Dirac points (K and K') are no longer independent of 
each other in armchair nanoribbon. Therefore, the num- 
ber of Dirac singularities is reduced to one in this system. 
On the other hand, the two Dirac points (K and K') 
are independent of each other in the case of nanotubc 
because the symmetric and antisymmetric combinations 
can satisfy the periodic boundary condition of nanotube. 

(VI) The phase space of armchair nanoribbon is given 
by 9 e (0,tt) and k G [-n/b,ir/b) [See Fig. 2(a)]. By 
comparing the phase space of armchair nanoribbon with 
the usual hexagonal Brillouin zone (BZ) of graphene, it 
can be seen that the former covers only the half region 
of the latter. To explore this feature further, we consider 
the other branch w- . Repeating analysis similar to that 
above, we obtain the corresponding energy dispersion re- 
lation and the wave function as 



e 2 = 1 + 4 cos 2 0' -4cos0'cos(A;Z), 



(21) 
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(a) 



fct 



2tt/6 



Phase Space Alternative BZ 1 




Alternative BZ 2 



Fig. 2. (a) The region denoted by shadow, 9 g (0, n) and k 6 
[— 7r/fe, tt/6), represents the phase space of armchair nanoribbon. 
The empty circle indicates the position of the Dirac singular- 
ity, (b) Two alternative Brillouin zones (BZ 1 and BZ 2) for 
graphene. The BZ 1 is useful to understand that the phase space 
of armchair nanoribbon is given by folding BZ 1 into the positive 
region of with respect to the fc-axis. The shadow region (BZ 2) 
is used to discuss on the Riemann sheets (see text). Note that 
alternative BZs contain two singularities corresponding to the K 
and K' points, while the phase space of armchair nanoribbon has 
a single singularity. 



b J _ a .-Q/-i) sin(Jfl') 1 
sin 9' \ 



2 cos6' 



— e 



where 9' e (— tt, 0) is defined by 



2 cos 19' 



w 



+ w- + \/z + z* - 2 + 4e 2 



(22) 



(23) 



By subtracting eq. (23) from cq. (12), we get cos 9' = 
cos 9 + cos(kl). Using this relationship between 9 and 9', 
we see that the point (k, 9) = (0, 27r/3) corresponds to 
9' = — 7t/3, for which cf) J _ reproduces (j) J + . Generally, by 
setting 6' = 6 — 7r in eq. (22), we obtain the wave func- 
tion 4> J + of cq. (16). Thus, the w_ branch does not yield a 
state which is independent of the w+ branch. Hence, we 
may consider the w+ branch only. This observation al- 
lows us to interpret the phase space of armchair nanorib- 
bon shown in Fig. 2(a) as the region given by holding the 
alternative BZ 1 of nanotube shown in Fig. 2(b) into the 
positive region of 9. In the remainder of the paper, we 
abbreviate ■ 



bi(s,k,9) as ^(k,9), 



ikl(J-l) 



■n sin J 9 ( e 



sin# 



-ie(fe,6») 



— s 



(24) 



It is interesting to note that the w + and cj_ branches 
correspond to different states in the case of nanotube. 
To show this, we consider two Riemann sheets for the 
variable z — e tkb . Thus, we consider the region, k € 
(— 2-rr/b, 2w/b], or the alternative BZ 2 in Fig. 2(b), and 
choose the branch w = \/z as a single valued function 
for the two Riemann sheets. From the energy dispersion 
relation eq. (13) we see that two points in the Riemann 
sheets, (0, 2tt/3) and (2w/b,n/3), result in e = 0. Note 
that k = and k — 2-rr/b belong to the different Riemann 
sheets. Hence, these points correspond to the indepen- 
dent Dirac singularities at the corners of the first BZ of 
a nanotube. 



3. Applications 

In §3.1, we examine recent experimental result of scan- 
ning tunneling microscopy (STM) by using the analytical 
solution eq. (24). 

3.1 STM Topography and Dirac Singularity 

The spatial dependence of the wave function (j) J s {k,8) 



is given by 



-iklj 



sin( J6*). Hence, the probability den- 



sity \4> J S (k, 9)\ 2 is proportional to sin 2 (J6*) and exhibits a 
spatial oscillation. Note that the probability density does 
not depend on the wave vector k. As a result, a STM im- 
age should be stable against a change in the wave vector 
k caused by, for example, a change of the bias voltage. 

As we have seen, the Dirac singularity corresponds to 
the state with 9 — 2ir/3 and k = 0. Therefore, the prob- 
ability density for a state near the Dirac singularity with 
9 = 2ir/3 disappears at J = 3i, where i is integer, as 
shown in Fig. 1(b). This behavior of the periodic absence 
of probability density (line nodes) has been recognized in 
literature. 8, 13, 16, 19 > On the other hand, for a state satis- 
fying 9 = 2tt/3 + 59, that is, for a state away from the 
Dirac singularity, the densities at J = 3i behave accord- 
ing to sin 2 (3i<5#) . Hence, as i increases or as the distance 
from armchair edge increases, the density at J = 3i in- 
creases gradually. 13 ^ Consequently, the behavior of the 
periodic absence of the probability density at J = 3i is 
not seen for the case that i is sufficiently large or away 
from the armchair edge. 

In the STM topography obtained by Yang et al. [Fig. 3 
of Rcf. 10], the density at J = 3i is suppressed for the 
case i < 3 and the appreciable density appears for the 
case i > 3. This feature of the STM topography allows us 
to interpret the character of their sample in the following 
way: the state observed by their STM was a state away 
from the Dirac singularity, that is, a state with 59 =/= 
0. If we assume that the probability density at J = 3i 
with i — 4 is comparable to that at J — 1, we have the 
following equality determining 59, 



sin 2 {1259) = sin 2 (27r/3). 



(25) 



The minimum value satisfying this equation is 59 = 
7r/36. This value of 59 corresponds to the states with 
energy \E\ « 7v / 3|c>6'| = 0.45 eV. Note that the energy 
value coincides with the energy difference between the 
Dirac point and the Fermi energy for graphene on 6H- 
SiC(0001) as pointed out in Ref. 10. Note also that a 
small 59 gives rise to a long-range electron-density oscil- 
lation of length (a/2)(ir/58)/3. A long-range oscillation 
with a period of about 2.5 nm is actually observed by 
Yang et al. The value of 59 which reproduces 2.5 nm is 
about 7r/60. The actual period of a long-range oscilla- 
tion may depend on the details of the substrate and the 
perturbations at the boundary. 

Here, we examine the effects of boundary perturba- 
tions, such as a change in bond length or a change in 
on-site potential, on the electronic state. First, suppose 
that the hopping integral between the A and B-atoms at 
J = 1 differs from that at J ^ 1 by 5-f. Then the first 
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order energy shift is given from cq. (24) by 

Ae = [#]t (5 Wx ) 4>l = -2s8-/\C\ 2 cos 9. (26) 

Because the energy shift is proportional to s, a high- 
est energy state in the valence band and a lowest en- 
ergy state in the conduction band undergo the opposite 
energy shift. This means that the energy gap changes. 
For example, an armchair nanoribbon with N = 3i — 1 
acquires a finite energy gap due to a change in bond 
length at the boundary. 13 ) Since Ae oc cos0, the gap is 
actually produced as a consequence of the shift, 66, of 
the Dirac singularity point along the axis of 9. Fujita et 
a?. 20 -* showed by a numerical calculation that a nonzero 
6j, that is, a nonzero 69 can be induced by Peierls in- 
stability for the case that TV is very small like 3 or 6. 
Since TV is large in the experiment, the contribution of 
the Peierls instability to 69 may be ignored. 

Next, we consider that the on-site potential at J = 1 
differs from that at J ^ 1 by v. The first order energy 
shift due to the boundary potential is given by 

Ae = (va ) <j>\ = 2v\C\ 2 . (27) 

Because the energy shift is independent of s, a highest en- 
ergy state in the valence band and a lowest energy state 
in the conduction band undergo the same energy shift. 
This means that a boundary potential does not change 
the energy gap. However, this fact does not necessarily 
mean that the boundary potential does not produce an 
observable effect in STM topography. The wave function 
can be sensitive to the presence of the boundary poten- 
tial. In fact, in the presence of v, the recurrence equation 
eq. (2) is modified as cf) 2 = —T~ X {K + vao)^ 1 , so that 
we obtain 

<f> 2 + = -(\++vT- 1 )cj> 1 + , (28) 

instead of eq. (15). By putting this new boundary condi- 
tion into eq. (14), we obtain the modified wave function, 



w. 



— w_ 



■ (J _i) sinJfl ! 

sin0 V+ 
(,/-2) sin(J - 1)9 
sin# 



vT 



(29) 



When v is sufficiently large in this equation, only the 
last term is important. By using the definition of T^ 1 in 
eq. (5) the last term can be rewritten as 



-w. 



-(j-i) sin(J- 1) 



sm( 



w + 
w + 



(30) 



It vanishes at J = 1 for any 9, which shows a change 
of the boundary condition. Note that the last term van- 
ishes at J = 3i + 1 for a state near the Dirac singularity 
(9 = 2n/3), while the original wave function vanishes at 
J = 3i for the state. It is also interesting to note that 
the boundary potential can change the direction of pseu- 
dospin. 

It is difficult to predict the value of boundary potential 
v, because there are a number of factors that can change 
the value of v. For example, v depends on functional 
group attached to the carbon atoms at the boundary. 
However, some information about v can be derived from 



an experimental STM image. To this end, we consider 
the amplitude at J = 2. Then, for a state near the Dirac 
point (9 = 2n/3), we obtain from eq. (29) that 

^.«^+«T-V+- (31) 

If \v\ — 0(1), a strong interference between the first and 
second terms can be expected in general. There is a possi- 
bility that the interference effect results in a strong differ- 
ence between the electron-densities at J = 1 and J = 2. 
Since the STM image of Yang et al. [Fig. 3 of Rcf. 10] 
shows clearly that the electron density at J = 2 is sim- 
ilar to that at J = 1, which indicates that \v\ is much 
smaller than 0(1). On the other hand, the existence of a 
small boundary potential is reasonable because a bound- 
ary potential about \v\ w 0.1 can be induced by intrin- 
sic perturbations, such as next nearest-neighbor hopping 
integral 21 ' and phonon. Detecting in STM images the 
signal of boundary potential v is an interesting subject. 

3.2 Raman G Band 

The analytical wave function is useful in calculating 
the clcctron-phonon (el-ph) matrix element for a process 
observed in Raman spectroscopy. Let us consider the G 
band in Raman spectra. It is known that two optical 
phonon modes at the Y point, whose displacement vec- 
tors are denoted by u x and u y shown in Fig. 3, contribute 
to the G band. The displacement vector u y is parallel to 
the armchair edge and u x is perpendicular to it. Here, 
we calculate the el-ph matrix elements for the u x and u y 
modes. 





Fig. 3. The displacement vectors, u x and u y , of the F point op- 
tical phonon modes. 67 (<57 ; ) represents a change of the nearest- 
neighbor hopping integral caused by the vibration. 



We define G x and G y as follows, 



G x 
67 



G 



N 



= ]T(^)t [Ta z ^ +a z T- 1 <f> J - 1 ] +c.c, 



J=i 

N 



j=i 



- {T(j) J+1 +T- 



1^-1 



} - <7x<i 



+ C.C. 



(32) 



It is straightforward to show that G x and G y are the 
expectation values of the perturbations caused by the 
lattice displacements u x and u y , respectively. Here, ^7 
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and Sj' are the changes of the hopping integral for the 
C-C bonds denoted in Fig. 3, produced by the lattice 
displacements u x and u y , respectively. By putting 4> J a of 
cq. (24) into (f> J of eq. (32), we obtain 



G x = 0, 



G 

5"/ 



V - = -3(cr x ) - e, 



(33) 



where (<j x ) is the expectation value of the Pauli matrix 
o x defined by 



N 

(a x ) = ^(</» s J )W s J = -scose. (34) 
J=i 

In eq. (33), we have obtained G y by putting the energy 
eigen equation of eq. (6) with K = ecro + cr x into the 
right-hand side of G v in eq. (32). 

We note that this result of G v derived by using tight- 
binding model is slightly different from the result of an 
effective-mass theory. The effective-mass theory gives 7 ) 

^ = -3K>- (35) 

The difference between the two schemes is e. Since e is 
scaled by the hopping integral, e vanishes in the con- 
tinuum limit, so that this e term does not appear in the 
effective mass theory. Here, we consider two direct conse- 
quences of the e term. First, because the Raman intensity 
is proportional to Gy, the e term may give rise to a small 
dependence of the Raman intensity on the laser energy. 
The dependence about 10%/eV is acceptable for the pro- 
cess that satisfies |(<r x )| w 1. Second, the strength of G y 
decreases for the state near the top (bottom) at the con- 
duction (valence) band. Thus, the Raman intensity for, 
if any, such a high energy state should be suppressed 
strongly. 

The result in eq. (33) means that the u x mode is not 
Raman active mode near the armchair edge and only the 
u y mode can be Raman active. This is a phenomenon 
peculiar to armchair edge 7,14 ) because it can be shown 
that the G x and G y do not vanish in the case of a periodic 
graphene as 

s (cos(e + e) - cos(e -e + kb)), 

(36) 

-3(<r x ) - e. 



c>7 



Recent experiments for the G band Raman spectroscopy 
for graphene edge by Cong et al. 22 ^ and Begliarbekov et 
al. 23 ^ support this prediction. 

4. Discussion and Summary 

It is often said that a step edge located on top of a 
highly oriented pyrolytic graphite (HOPG) can have a 
regular structure over a large length more than 10 nm, 24 ' 
The wave function for such a step edge can also be con- 
structed analytically. The equation for a step armchair 
edge on top of HOPG is given by replacing the matrix 
K in eq. (6) with K — ma z . Here, m represents the po- 
tential difference between the A and B sublattice which 



is induced by the stack in the Bernal configuration. It 
is easy to show that the spinor eigenstate of cq. (9) is 
modified by the presence of m as 



'l-z + V(-l+^) 2 +4z(e 2 



2 ) 



2z{e—m) 
1 



(37) 



Note that the lowest energy state with k = (z = 1) in 
the conduction band has amplitude only on the A-atoms 
for the case that m > 0. It is obvious that by taking the 
limit 6 — )■ tyl in eq. (37), the normalized spinor is written 

as 



<t>+ = 



(38) 



which is a pseudospin polarized state. The electron- 
density pattern of this state is shown in Fig. 4, where 
the density forms the hexagonal pattern. It is important 
to recognize that the appearance of the hexagonal pat- 
tern results from two causes: (1) the asymmetry of the 
potential energies at the A and B sublattice and (2) the 
existence of armchair edge. The condition (1) is satis- 
fied in the bulk of HOPG 25 ) and epitaxial graphene on 
SiC(OOOl). 26 ) The STM images show a triangular lattice 
there because the amplitude can be observed only at one 
of the two sublattice for the case that m ^ 0. In addi- 
tion to (1), near the armchair edge (2), the low energy 
electron does not have amplitude at the sites J = 3i [sec 
Fig. 1(b)] due to sin(J0). As a result, a low bias STM 
topography observes the hexagonal electron-density pat- 
tern near a step armchair edge on top of a HOPG. 11 ' 27 ) 




Fig. 4. The electron-density pattern of a low-energy state near 
the armchair step edge located on top of a HOPG. The unit of 
(\/Za X \/Za) 30° superstructure is shown by the rhombus. This 
unit cell is the same as that of the Kckule pattern. 



Beside the hexagonal pattern, an electron-density pat- 
tern which is referred to as (v^a x y/3a) 30° or (\/3 x 
\/3)i? 30° superstructure is frequently observed near de- 
fects on graphite, such as metal particles, 28 ' 29 ) isolated 
adsorbed molecules, 30 ' 31 ) step edges, 24 ' 25 ' 32 ) and lattice 
vacancies. 25 ' 33 ~ 35 ) This superstructure is denoted by the 
rhombus in Fig. 4 (see ref. 28). The periodicity of the 
(y/3a x y/3a) 30° superstructure is actually included in 
the wave function for armchair edge as shown in Fig. 4. 
Note, however, that the origin of this superstructure 
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observed in the experiments for step edges can not be 
attributed only to armchair edge since the presence of 
zigzag edge might play an important role. A simulation 
by Niimi et al. 27 ^ indicates such a possibility. To fully 
specify the origin of the (y/3a x y/3a) 30° superstruc- 
ture, it is desirable to construct analytical wave function 
for a mixed edge consisting of zigzag and armchair edge 
parts. 

In summary, we have constructed the analytical so- 
lution of the tight-binding model for armchair nanorib- 
bon. The wave function is described by the superposi- 
tion of the incident wave to the armchair edge and the 
scattered wave, i.e., the formation of standing wave. The 
scattering process is intervalley scattering, in which the 
pseudospin is invariant. The difference between armchair 
nanoribbon and nanotube appears with respect to the 
number of Dirac singularities in their phase spaces. Since 
the armchair boundary condition identifies the half of 
graphene BZ (including the K point) and the other half 
of the graphene BZ (including the K' point), the phase 
space of armchair nanoribbon is given by the half region 
of the graphene BZ. As a result, in the phase space of 
armchair nanoribbon, there is a single Dirac singular- 
ity whose wave function is written by the antisymmetric 
combination of the two states at the K and K' points. We 
have examined the recent STM topography by using the 
analytic solution. It seems that the STM image obtained 
by Yang et al. 10 ^ corresponds to a state away from the 
Dirac singularity. This speculation is consistent with the 
appearance of a long-range oscillation in their STM im- 
age. Moreover, the observed STM image does not allow 
assuming the existence of a strong boundary potential 
of order of the hopping integral. A lattice distortion at 
armchair edge can be taken into account as a shift of the 
Dirac singularity point along the #-axis, which is similar 
to the curvature effect in metallic carbon nanotubes. In 
addition, the analytical solution allows us to observe eas- 
ily that u x mode is not Raman active near the armchair 
edge. The analytical solution simplifies greatly the anal- 
ysis which was done numerically with the tight-binding 
model, 14 ) and is also useful for confirming the result ob- 
tained analytically with the effective-mass model. 7 ) 
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